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Abstract 

We study distances of propagation and the group velocities of the muon neutrinos in the presence 
of mixing and oscillations assuming that Lorentz invariance holds. Oscillations lead to distortion 
of the f^— wave packet which, in turn, changes the group velocity and the distance z/^ travels. 
We find that the change of the distance, dose, is proportional to the length of the wave packet, 
' (Jx, and the oscillation phase, (pp, acquired by neutrinos in the it— and K— meson decay tunnel 

where neutrino wave packet is formed: dose oc a (pp. Although the distance dose may effectively 
correspond to the superluminal motion, the OPERA result can not be explained. For OPERA we 
r~| . obtain dose ~ — 10~^ cm using Am^ = 2.5 x 10^^ eV^. We analyze various possibilities to increase 

dose and discuss experimental setups in which dose (corresponding to the superluminal motion) can 
reach an observable value ~ 1 m. 
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I. INTRODUCTION 



The OPERA collaboration has measured velocity of the muon neutrino which apparently 
exceeds the velocity of light in vacuum Earlier MINOS gave some indications of su- 
perluminal motion of neutrinos though with lower statistical significance 0]. Plethora of 
different interpretations and possible explanations of these results has been proposed . A 
series of papers Q, [sj, joj, jsl employed an analogy with optics where the group velocity 
larger than the velocity of light was predicted without Lorentz violation and later observed 
in several experiments (see e.g., [HI, [12I). The effect is related to distortion of 

the photon pulse during propagation in media, so that suppression of trailing edge of the 
pulse leads to increase of the group velocity. This superluminal motion does not contradict 
causality and no information can be transmitted with velocity larger than the velocity of 
light. 

In 043] the superluminal motion of the muon neutrinos has been considered in the 
presence of flavor mixing, when muon neutrino is described by coherent combination of the 
mass eigenstates. Distortion of the z/^ wave packet is produced by interplay of two effects: 
(i) the coordinate independent z/^ oscillations in time, and (ii) a relative shift of the wave 
packets of mass eigenstates due to difference of their group velocities ^. In this proposal 
an additional contribution to the distance of t'^ propagation, d, and to the group velocity 
of muon neutrinos is proportional to the difference of the group velocities of the neutrino 
mass eigenstates Av ~ Am'^ /2E^, which is extremely small. Since d ~ the effect can 

be enhanced if the survival probability of muon neutrino, P^^, is very small. The latter, in 
turn, requires mixing to be very close to the maximal one and the oscillation phase to be 
very close to vr. Apparently, this can not explain the OPERA result for several reasons: (i) 
Small survival probability means that only small number of events should be observed in a 
detector, in contrast to observations, (ii) Fine tuning of the phase can be achieved at certain 
energy of neutrinos only, and introduction of new matter effect [l3| does not resolve this 
problem, (iii) Additional distance of Vjj, propagation is restricted by the length of neutrino 
wave packet which is much smaller than 20 m indicated by OPERA. 

In this paper, we study the group velocities of neutrinos in vacuum and in matter within 
the framework of Lorentz invariance. We show that the dominant effect, which modifies 
velocity of the muon neutrino, is a distortion of the wave packet by the coordinate- 
dependent oscillations. That is, by the oscillations within the size of z/^ wave packet in the 
coordinate space. The effect of relative shift of the wave packets discussed in the literature 
is negligible. We show that the oscillation distortion of the wave packet still can not explain 
the OPERA result: an increase of the distance of propagation, dose-, is much smaller than 
the required 20 m. Furthermore, dose rapidly increases with decrease of the neutrino energy 
in contrast to observations. 

We study in detail properties of additional contributions to the distance traveled by neu- 
trinos and to the group velocity. We estimate the distance d for various existing experiments 
and indicate new experimental setups where the effect can be observed. 

The paper is organized as follows. In Sec. II we present general formulas for the distance 
traveled by oscillating neutrinos during a given interval of time. In Sec. Ill we construct a 
wave packet of the muon neutrino produced in the pion decay and describe its properties. 
In Sec. IV we give explicit expressions for the distances traveled by muon neutrinos. We 



We will give detailed explanation of the effect in the Appendix. 
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consider the limit of small size of decay tunnel when the wave packet has approximately the 
box-like shape. Effect of relative shift of the wave packets of mass eigenstates on the distance 
is computed. In Sec. V we estimate the distances and velocities for existing experiments and 
propose experimental setups in which the effective superluminal motion might be observed. 
In Sec. VI we consider group velocities of neutrinos in matter. We conclude in Sec. VII. In 
the Appendix a simple explanation of the shift effect proposed before is given. 



II. OSCILLATION PROBABILITY AND DISTANCE OF i/^ PROPAGATION 



For simplicity we will consider the two neutrino mixing 



v. 



CZ/i + SZ/2, I^r = CZ/2 - SZ/i, (l) 



where c = cos 6', s = sin 6, and 6 is the mixing angle; i>i {i = 1,2) are the mass eigenstates 
which correspond to eigenvalues m,. We will show in Sec. |Vl] that this is a good approx- 
imation for the three neutrino mixing case. Recall that in OPERA, MINOS, T2K, etc., 
experiments the neutrinos oscillate in matter and therefore in the presence of non-zero 1-3 
mixing one should use the mixing angles and mass splitting in matter. However, apart from 
the resonance regions E ~ 0.1 GeV and E = {A — 8) GeV the 2z/ description with vacuum 
values of mixing angles and Am^ = — ml gives very good approximation. 

Evolution of the muon neutrino state after it exits the region of formation of wave packets, 
i.e. a decay tunnel, is described as 

\u^{t)) = cfi{x - vit)e~'^''+'P'''\ui) + sf2ix - V2t)e-'^''+'P''''\u2), (2) 

where are the shape factors, Ei and pi are the average energies and momenta of the wave 
packets, respectively, and Vi are the group velocities of the mass eigenstates. The shape 
factors are normalized as 



/ 



dx\Mx - v,t)\^ = 1. (3) 



If muon neutrino is detected at some moment of time t, the wave function of z/^ (the 
amplitude of probability to find z/^) is given by 

iP.^{x,t) = (z/^|i/^(t)) = cVi(x-t;it)e-^^i*+*P^^ + s2/2(x-t;2t)e-^^^*+'P^^ (4) 

Here we used implicitly that the shape factor of detected muon neutrino is constant in 
the detection area and does not depend on time. Alternatively, one can consider fi as the 
effective shape factors which already include process of detection. 
The survival probability at a moment of time t equals 

P,,{t)^ j dx\ij,^{x,t)\\ (5) 

In what follows, we will compute the averaged coordinate {x{t)) of the muon neutrino in 
a moment of time t defined as 

= Idxx\^,^{x,t)\' 
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where the denominator is nothing but P^^{t). 
To a very good approximation we take 

/iW = /2(x) (7) 

at the initial moment of time so that difference between /i and /2 in an arbitrary moment 
of time arises solely due to the difference of group velocities of vi and V2- The difference 
/S.V = V1—V2 produces a relative shift and eventually separation of the wave packets of mass 
eigenstates in the process of propagation. The separation leads to loss of coherence between 
the mass eigenstates. Numerically, 

Av = vi-V2^ —— = 5 X 10"^^ ( z — K II- 



Using the expression (jl]) and the normalization condition ([3]) we find 

P,,{t) = c' + s^ + 2ch'l, (8) 

where 

/ = y dx fi{x — Vit) f2{x — V2t) cos^{x,t), (9) 
and the relative (oscillation) phase of mass eigenstates, $(x,t), is given by 
^{x, t) = AEt - Apx, AE = El- E2, Ap = pi - pa- 
The phase $(x,t) can be split into two pieces as 

$(x,t) = + X, (10) 

where 

(f) = AEt - Apvt, V = (11) 

is the standard oscillation phase which depends on time only and 

X = -Ap{x -vt). (12) 

The phase x depends also on the coordinate x and, as we will see, describes the change of 
the total oscillation phase within the wave packets. 

Let us introduce the average coordinates of the mass eigenstates as 

(xi) = j dx X [fi{x - Vit)]'^ =Vit + j dz z [fi{z)f. (13) 

If fi{z) is symmetric with respect to the central point of packet, Zq. fi{zQ — h) = /{{zq + h), 
we obtain (performing shift of integration z ^ z — Zq): {xi) = Vtt + Zq. 

Using (xi) we define the averaged coordinate of the two mass eigenstates, 

1 If 

^ = ^ii^i) + {^2)) = 2(^1 + V2)t + dz z [fiiz)]\ (14) 
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and the half-difference of the coordinates 



Ax 



2((^i) 



(X2)) = ^Avt. 



Then, the numerator of {x{t)) in ([B] 

dx X \ipunix,t)\'^ 



can be written as 
C + S )x+ - 



+ 25^0 V, 



where 



J 



dx X fi{x — Vit)f2{x — V2t) cos$(a;,t) 



(15) 
(16) 



is the average coordinate of the overlap of the two mass eigenstates. 
Using definition (Q, ([8]) and (fT5|) we can present {x(t)) in the form 

Avt cos 29 sin^ 29 



{x{t)) 



X 



-{J-xI). 



(17) 



171) is the distance traveled by massive neutrinos, and the two others 



2P 2P 

tin. I 

The first term in 

describe the additional distances due to oscillations and relative shift of the wave packets. 
The distance of z/^ propagation during a time interval from to to t equals 

d{t) = {x{t)) - (x(to)). (18) 

As to we will take the moment of time when the neutrino wave packets are completely 
formed. For neutrinos from pion decay to = Ip/ v-j^ is the moment of time when pion reaches 
the end of decay tunnel. According to f|T71) the distance d(t) can be presented as 

dit, to) = dlight + dmass + dosc-a + dghift-s, (19) 

where 

dlight = c(t — to); 

is the distance traveled by the light. 



d„ 



m 



mi + mi . 







(where E is the neutrino energy) is the contribution to the distance due to non-zero neutrino 



mass; 



dn 



sin^ 29 



J-xI 



P. 



J-xI 



P. 



till 



to 



is, as we will see, the contribution to the distance from oscillation distortion of the muon 
neutrino wave packet and from shift of the packets in the case of asymmetric shape factors; 
and 



d 



Av cos 29 



shift— s 



t 



to 



LP,,it) 



(20) 



is the contribution from the relative shift of the wave packets and coordinate-independent 
oscillations in the case of symmetric shape factors. Notice that the terms in the brackets 
evaluated at t = to are negligible, since usually t ^ to and P^^ito) ~ 1. However, it may 
not be always the case for dosc-a because (J — xl) is non-linear function of t. 

As we will see, the distance traveled by muon neutrinos is an oscillatory function of time. 
Correspondingly, the velocity changes with time. Therefore we can speak about the average 
velocity for a given time interval {t — to) defined as v{t) = d{t, to)/{t — to). 
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III. NEUTRINO WAVE PACKETS FROM PION DECAY 



For definiteness we will consider neutrinos from pion decays. Actually, these neutrinos 
dominate in the neutrino fluxes of MINOS, T2K and OPERA experiments. The contribution 
from K— mesons can be considered similarly. We describe here the wave packets of muon 
neutrinos (the wave function -0, , , , in the eq. (jlj)). The strict derivation of the expressions 
for the wave packets is given in [m] , and here we explain properties of the shape factors and 
phases using simple physics arguments. 

Let us consider first the shape factors. Pions are produced in the strong interaction of 
protons with a solid state target. Protons have wave packets of very small size and can be 
considered as point-like. The target nuclei are localized within the atomic scale. Therefore 
we can take in our computations that pions are produced at fixed space-time point x = t = 0. 
The wave packets of pions are very short in the configuration space, so that pions also can 
be considered as point-like with space-time trajectories of motion x = f^t. Here is the 
group velocity of pion. In experiments under consideration pions are ultra-relativistic with 
Lorentz factor 7^ = i?^/m^ 3> 1, where E.,^ and are the pion energy and the mass, 
respectively. 

Pions decay in a decay tunnel of size Ip which includes the length of the horn area and 
the decay pipe. Essentially Ip is the distance from the target to the end point of decay 
pipe where survived pions are absorbed. We neglect the effects of interactions of pions with 
particles in the decay pipe. In the decay tunnel the quantum states are a superposition 
of un-decayed pion state and the state of muon and muon-neutrino. The neutrino waves 
emitted from different points of the pion trajectory are coherent, thus forming a single wave 
packet. 

The decay length of pion equals 

where and r° are the lifetime of pion in the laboratory frame and the one in the pion rest 
frame, respectively. r° = Pq ^ with Pq being the pion decay rate at rest. 

If Ip < Idecay (high energies, long-lived parent particles, short decay tunnel), the wave 
packets of mass eigenstate z/j have the size 

when neutrino is emitted in the forward direction with respect to the pion velocity. The 
wave packet size is proportional to Ip] the factor in bracket represent shrinking of the packet 
due to pion motion jl5|. In the second equality in fl^T]) we used <^ m^r. Numerically ^ 

a ~ 50 cm f^V' ( = 17.7 cm ( ' ( . (22) 

If Idecay <^ (low energies, long decay tunnel, short-lived parent particles), the size of 
the wave packet is determined by the decay length: 



Idecay ^ " ^^)- (23) 
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For two body pion decay the neutrino energy E is determined by Lorentz factor of pion, , or vice versa: 

E ^ 
IGcVV 



77r — 16.8 ( .^y ) to a good approximation for E > 500 MeV and forward going neutrinos. 
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For relativistic pions (7^ > 1), d ^ Ir^-f-^ instead of In the non-relativistic limit, 

Vt, — )■ 0, the size becomes the largest: cij = r°t>j ~ 7.8 m. 
Thus, in general, 

CTj = lform I ] ^ , (^4j 

/ 272 

where //orm is the region of formation of the neutrino wave packet: 

lform ~ min {Ip, I decay} ■ 

In what follows we will neglect difference between the sizes of wave packets of the mass 
eigenstates, taking ai ~ (J2. It is an excellent approximation because the difference is of the 
order ~ IpAv ~ IQ-^'^cm for Ip = 1 km, E = 1 GeV and Am^ = 2.5 • IQ-^ eV^. 



The neutrino shape factor from pion decay can be written as [14 



/(a;) = /oe^(^-'^)n(x,[0,a]), (25) 

where 11 [x, [0, o"]) = 1 in the interval < a; < a and 11 vanishes outside this interval. The 
box function reflects a finite size of neutrino production region with sharp edges related 
to the definite points of pion production and absorption at the end of decay tunnel. The 
exponential factor follows from the amplitude of the pion decay with the rate enhanced by 
the factor l/(f — f^r) which corresponds to shrinking of the neutrino wave packet emitted in 
forward direction (see eqs. This is nothing but the Doppler effect. 

The normalization condition ([3]) gives 

where 

y = ^^ = ^. (27) 

'^■K I' decay 

If Idecay > Ip, wc have lform = Ip, SO that y is the size of the decay tunnel in the units of the 
decay length. Numerically, 

^ ^P^O _ 1 00 /^77r\"W Ip 



1/ = -^^:^^1.28x ^ {-^]. (28) 

Idecay 7. VIO/ VlOOm; ^ ^ 

If Idecay < Ip and wc do not consider effects in the exponentially suppressed tails, then 

lform ~ Idecay, SO that |/ ~ 1. 

In the approximation Vi ~ v, the difference of momenta of the neutrino mass eigenstates 
equals fl3 | 

Am2 Am^ , , , 

The factor in the denominator of (1291) is the same as the one which describes shrinking of 
the wave packet size in comparison with size of the wave packet formation region. The 
shrinking is accompanied by increase of frequencies by the same factor. This is again the 
Doppler effect related to neutrino emission from moving pion. Consequently, the oscillation 
effect within the packet is enhanced: It is determined by the oscillations within the wave 
packet formation region, rather than within the size of the packet itself in the laboratory 
frame. 
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IV. DISTANCES OF PROPAGATION 

We will use the formulas derived in Sec. [11] and the wave function of the muon neutrino 
constructed in Sec. [1111 to find P^^ and {x(t))osc- We first compute the effect of oscillations, 
which gives the dominant effect, neglecting a relative shift of the ui and 1^2 wave packets 
due to difference of group velocities. Then, we will estimate the effect of the shift neglecting 
oscillations along the wave packets. 



A. Oscillation effect 

Let us set vi = V2 = v {Av = 0) and compute I and J given in ([9]) and (fT6|) . We use the 
shape factors (!25|) and the phase $ with defined in (ITT!) and x given according to ( !T2|) and 
(UnDby 

Am^ X — vt 
^ ~ 2E v-v^' 

We obtain the probability 

P^^{t) = + + 2s^c^ {cos (f)Ic - sin0/,), (30) 

where 

/y 1 
dzf{z) cos\ Ap\z = , ,0 {-ye^y + y cos (f)p + sin (pp) , 

/y 1 

dzf{z) sin \Ap\z= ^ _ _^ (0pe"'^ + ?/sin0p - 0pCos0p) . (31) 

Here 

Am^ (T Am^ 

'^^ " ^ (^7^ " = 

is the phase change within the wave packet. Here li, = AtcE/Aitl'^ is the oscillation length. 
The second and third equalities in eq. fl52]) are valid for an ultrarelativistic pion: ^ 1. 
Thus, (pp is given by the phase acquired within the distance of formation of the neutrino 
wave packet, I form- For Ip < I decay we have 

Am2 , , / Am2 \ f E YW L \ 

— 2.5X10- (^^)(—) (^). (33) 

For Ip > Idecay inserting a from ( l23l) into ( l32l) we obtain 0p = 2t['~^t,t^v /l^. 

We now calculate (x(t)) ffTTj) for Av = 0. The average distance traveled by mass eigen- 
states f[T^ with the wave packets ([25]) equals 

a; = t;t + af-^-iV (34) 



1-e-y y 

Straightforward computations of the contribution to the distance of Up propagation from 
the oscillation effect give 

{x)osc = - 9^,9 (f^s sin (p- He cos (p) . (35) 

2Ppf, y^ + (p; 
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Here the coefficients Ks and Kc in front of sin and cos equal 

y 



1 - e-y 



-2y(j)pe ^ + g sin 0p — n cos 0p 



1/2 + 02 



1 



1- e-y y 



bpC y + y sin 0p - 0p cos 0p) 



(36) 



y 



(?/2 — 02)e y -\- q cos 0p + n sin < 



y' + 



1 Iw ^ 

1 — e y y / 



+ ?/ cos 0p + 



, sm I 



(37) 



where 



q = y\y-l) + (j)l{y + l), 



n 



My' - % + 0p) 



In reahstic setups: ?/ ~ 1 and 0p ^ y. For 0p — )■ we obtain (x)osc = 0, because there is 
no room in which oscillation effect develops if Ip ly. In the limit of small 0p the distance 
equals 



sin^ 29 ac 



2P 



y' 



Tjs sm ( 



Tjc — COS ( 

y 



(38) 



where 



1 - 



1 - e-y 



2 



For y = I the value ?7s/z/^ ~ 1/12.71 and it depends very weakly on y. It equals exactly 
1/12 in the case of box-like packet which corresponds to ?/ — )■ (see Sect. IVB). 

The effect of oscillations on the group velocity of the muon neutrino is illustrated in the 
Fig. [H In the case b), the oscillations suppress the front edge of the z/^ packet stronger than 
the trailing edge. As a result the "center of mass" of the packet shifts toward backward 
and therefore the distance of propagation reduces. In the case c) realized in other moment 
of time the trailing edge is suppressed stronger, the center of mass shifts forward and the 
distance of propagation increases. 



B. Limit of small production region, y 1 

Expression for (x) obtained in Sec. IVA substantially simplifies in the limit y —> 0. It 
corresponds to short decay tunnels Ip ^ Idee, or to high energy parts of the neutrino spectra 
in real experiments. In this limit F — )■ and therefore 

f{x) = ^U{x,[0,a]), (39) 

and we will refer to this as to the "box-like" packet. 
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dosc< 




dose > 



FIG. 1: The wave packets of the muon neutrino from pion decay without oscillations (a) and with 
oscillations in two different moments of time (b),(c). The lower (green) parts of the shape factors 
show the fraction, whereas the upper (blue) parts correspond to the Vr fraction which appears 
due to oscillations. The arrows indicate positions of "centers of mass" of the parts. The red 
boxes show shifts of the centers due to oscillations with respect to the center in the no-oscillation 
case. The pannel (b) corresponds to the baselines < L < /osc/2, when the front edge of the wave 
packet is suppressed. The pannel (c) is for /osc/2 < L < lose, when the trailing edge is suppressed. 

To obtain P^^ and {x)osc for the wave packet (!39l) we take the limit ?/ — )■ in ( l30l) and 
while keeping 0p finite. The survival probability becomes 



p^^(t) = c^ + / + sin^2^ 



sm 



cos 



2 



(40) 



In the limit ?/ — )■ the second term in (IMl) tends to cr/2, and consequently, x 
Finally, {x)osc takes the form 



vt + 



2 • 



(x(t)). 



sin^ 26 a 



2P 



MM 



2 sm — 
2 



f)p cos — 



sm 



+ 



(41) 



If 0p is small (for existing setups 



10" 



{x)osc varies with distance traveled as ~ sini 



For short baselines, when < 7r/2, the additional distance {x{t))osc < 0, i.e. oscillations 
suppress velocity (see Fig. [1]). In addition, for small the distance is suppressed by 0. 
At around tt/2, 3tt/2, etc. the oscillation effect is maximal for {x{t))osc- Furthermore, 
at ~ 7r/2 (37r/2) the oscillation probability is a decreasing (increasing) function of x. 
Therefore, the effective shape of the z/^ wave packet deforms in such a way that the effective 
velocity is smaller (larger) than the normal velocity of massive particles at ~ 7r/2 (37r/2). 
In the limit 0p — ?■ 0, the expressions in (HOi) and (jUJ) give 



{x{t))osc) 



c^ + s^ + - sin^ 29 



sin^ 26 



a0„ sm 



cos 



+ 



2 



COS( 



(42) 
(43) 



where we have used that the factor inside the parenthesis in front of the sine factor in (jH 
is approximately equal to ]^0p. We have not expanded the cosine and sine factors in (Il2 
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(H3|) since or its deviations from nvr could become smaller than 0p at certain moments of 
time. 

According to the contribution to the distance is proportional to ^ 

a<P, = 1.25 X 10-cm ( ( " f ^) " f ' . (44) 

Thus, in the limit Av — ?■ the total distance traveled by neutrinos equals 
d{t, to) = v{t - to) + [{x(t))osc - {x(to))osc] , 

where the first term is the distance traveled by a massive neutrino. The second term takes 
the largest positive value ^ 2 x 10~^cm, at around ~ ^ and for the typical values of the 
parameters in the parentheses of eq. (jS]). This term dominates over the first one. 

C. Estimation of effect of shift of the wave packets 

For realistic experimental setups Avt <^ a < Idecay The difference of group velocities 
of the mass eigenstates produces a relative shift of their wave packets and therefore an 
additional distortion of the wave packet. Due to the shift, three different spatial parts of 
the z/^ wave packet appear (we assume that vi > V2): 

1) . The front edge part: {v2t + u) < x < (fit + a). In this part only the Vi packet is 
present, therefore ip^^ ~ c^/(o") = c^/o and \ipu^\'^ = c^/o- The length of this part equals to 
Avt and we can safely neglect the change of f{x) within this interval. 

2) . The overlapping part: (fit) < x < (f2t + a). Here both wave packets are non-zero 
and interfere between each other: 

3) . The trailing edge part: (f2^) x < (fit). Here only the 1^2 wave packet is present 

and 4)^^ ^ 5^/(0), so that \ip,yj'^ ~ s'^f'^e'^^ . 

Apparently, oscillations do not affect contributions from the front and trailing edge parts. 
These parts produce additional asymmetry of the v^j, wave packet which is proportional to 

(e' - /a-^) A„t = (c- - .'a-^) (^) ^. 

Consequently, the contribution to the distance of propagation dshift-a ~ ^vt. As we 
will show later the asymmetry produced by the overlapping region (which depends on os- 
cillations) is even smaller than Therefore in the first approximation we can neglect 
dependence of the oscillation effect inside the wave packet size on distance and take the 
phase $ to be constant. 



^ Though and E are related with each other in pion decay we represent the both factors independently 
for possible use of the formula for neutrinos from other sources such as muon decay, etc.. 



(? + S^COS$(x) ( 1 + 



Aftr 
2(f - f.) 



(45) 
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Let us introduce the dimensionless parameter 

AvtT Avt 1 



2t'7r 2 Idecay 



(46) 



which is the half-shift of the wave packets in the units of the decay length. Its value can be 
estimated as 







2 fyrTTT 



When the oscillation effect in the coordinate space within the wave packet size is neglected 
we have J = coscp J dx x fi{x — Vit)f2{x — V2t) and / = cos0 / dx fi{x — Vit)f2{x — V2t). 
Then straightforward computations give 



^ = cosv 

1 — e y 



J = cos 



1 - e-y 



y 



(48) 



For the average distance traveled by the mass eigenstates we have x given in ( l34l) with 
V = 0.5(fi + f2). 

Using (HHil we obtain correction to the distance of z/^ propagation 

, . ^^ , sin2 26' f Avt\ ^, , 

{X{t))shift~a = COS0 



where 



i^(2/, 



'1 - 6-2^)2 



ye-^^ - (1 - e-^)(e-^ + e^"^) 



(50) 



and the probability P^^ is given in ([8]) with / from fl48l) . 

As can be seen from fHTj) . e is negligibly small for typical values of experimental param- 
eters. In the limit of small e we have 

= (13^ [2ye-^ - 1 + e-^y] . (51) 

Also in most of the experimental settings ?/ ~ 1, and consequently F{y,e) ~ In 
particular, F{1) = 0.322. Then according to fHOl) \d{t)\ ^ {x{t))shift-a ~ At>t, assuming 
COS0 ~ 0(1) and P^^ ~ 0(1). Notice that {x(to)) term in (i(t) can be ignored for baselines 
L > 100 km, since Ip <^ L. Numerically, 

A.t.^L^0.5xlO-(^^j(^j (^)cm, (52) 



and L = ct stands for the baseline distance. Therefore, contribution to d{t) due to velocity 
difference of mass eigenstates is extremely small. 

For small y we obtain F ^ — y/3. In the limit of ?/ — )■ 0, which correspond to long lifetime 
compared to the travel time inside decay tunnel, the function F and therefore the distance 
vanish: {x{t))shift = 0. This limit corresponds to symmetric (box-like) wave packet. So, 
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the contribution to the distance (H9|) described here corresponds to non-symmetric shape 
factors. 

If ?/ — )■ oo, we obtain F — )■ —1. Let us underhne that F < 0, and therefore the correction 
is negative (reduces the effective velocity). It can be positive if formally F < 0, that is, when 
the amplitude of the shape factor is a decreasing function of x, or equivalently, increases 
with time. 

Let us now estimate corrections to the distance due to the overlapping part by taking into 
account the coordinate-dependent oscillation effect. The correction due to At>t in fHSj) does 
not depend on x, and therefore can be "absorbed" in redefinition of the mixing parameters: 

I V-., r ^ (1 + rs^ff{xf \cl + s\ cos ^{x)\\ 

where 



and 



1 1 + rs^ ' ^ 1 + rs^ ' 
/\vtT Avt 



2{v -v^) a ' 
According to fHOl) the resulting additional distance equals 

{x)sMSt-a oc + p'^'"'"' = (1 + r)s\\ 



Therefore the correction to distance of propagation due to the shift is 



cr 

where {x)^i^if^_a is the shift for constant survival probability. 

V. NEUTRINO VELOCITY IN VARIOUS EXPERIMENTAL SETTINGS 
A. Distances for existing experimental setups 

Here we present numerical estimates of the correction terms in d{t) ignoring the length 
of pion decay tunnel compared to baseline distance, so that d(t) ~ {x(t)). According to (fT9|) 
the corrections to normal distance traveled by massive neutrinos consists of three different 
terms. It may be instructive to compare behavior of these terms in the limit of small phases 
<^ 1 and 0p <^ 1 when all the contributions become linear in the time interval: 

sin^ 29 a(j)r) Am^ , , sin^ 2Q aE(j)„ ^ , 

dose = TT — it-tn) = — -^Av(t-to), 

2P^^{t) 12 2E ^ °^ 2P^^{t) 12 ^ 
cos2^ . . 
dshift-s = 2P^^{t) 

sin^ 26 

dshift-a = ^^-^F{y,e)Av{t-to). (53) 
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TABLE I: The values of y = Iform/ldecay, sine of the oscillation phase, (j), the wave packet length, a, 
as well as the contributions to the distance of i/^ propagation from the mass terms, dmass, from the 
relative shift of the wave packets, dghift-s, and from oscillations, dose- We use Am^ = 2.5x lO^^eV^, 
and sin2 29 = 0.97. 



Experiment 


E (GeV) 


y 


sin (f> 


a (cm) 


-dmass (cm) 


dshift-s (cm) 


dose (cm) 




OPERA 


17 


0.491 


0.270 


0.671 


1.58 X 10" 


-16 


2.80 X 10" 


-17 


-1.0 X 10" 


5 


OPERA 


1 


1 


-0.99 


23.3 


4.56 X 10" 


-14 


1.65 X 10" 


-14 


+1.6 X 10" 


3 


MINOS 


3 


1 


1.00 


7.7 


5.10 X 10" 


-15 


1.70 X 10" 


-15 


-4.6 X 10- 


4 


MINOS 


1 


1 


-0.995 


23.3 


4.59 X 10" 


-14 


1.62 X 10" 


-14 


+1M X 10- 


-3 


T2K 


0.6 


1 


0.052 


38.8 


5.12 X 10" 


-14 


2.94 X 10" 


-13 


-2.2 X 10" 


3 


T2K 


0.4 


1 


-0.997 


58 


1.14 X 10" 


-13 


4.25 X 10" 


-14 


+3.9 X 10" 


3 



It follows immediately from these equations that dose S> dn 



d 



shift- 



The contribution 



dose is strongly enhanced by factor aE in comparison to the other contributions. 
All the terms in f|T9|) quickly decrease with energy: 

1 . 1 . 1 



,_e.^. oci. 



The contribution dshift vanishes for maximal mixing. Other contributions rather weakly 
depend on (or independent of) 9. 

In the Table [T] we present numerical estimate of the correction terms d{t)i ^ {x)i for 
different existing experiments. We take Am^ = 2.5 x 10~^eV^, and sin^ 2^ = 0.97, which 
corresponds to the largest departure from the maximal mixing allowed at 90% CL by the 
SuperKamiokande atmospheric neutrino data [ig'] We use eqs. (135|) and (!38|) to compute 
dose- To find d^ass "we have assumed that m'^+m\ ^ Am? = 2.5 x 10~^eV^ which corresponds 
to the strong mass hierarchy, dghift-s has been computed according to eq. (|2U]) . Notice that 

dshift—s 

packet 



Mmass I cos 26'/P^^. As follows from (H9l) for the exponential (asymmetric) wave 

F{y) sin^ 20 cos0 



d 



shift— a 



d 



shift- 



2 COS 2^ 

and therefore dghift-a ~ dshift-s unless cos 2^ is very small. 

In the Table I we show also the size of the neutrino wave packet a which gives the absolute 
upper bound on additional contribution to the distance of propagation, the parameter y 
and the sine of oscillation phase. The phase 0p acquired by neutrinos in the wave packet 
formation region equals 0p = 4 x 10~^ for E = 17 GeV in OPERA experiment and 0p ~ 
8.2 X 10"^ in all other setups with y ~ 1. The contributions have been computed for 
the average energies and for representative low energies in the spectra in corresponding 
experiments. We used the following baselines and decay tunnel lengths: OPERA: L = 730 
km and Ip = 1095 m, MINOS: L = 735 km, and Ip = 715 m, T2K: L = 295 km and Ip = 118 
m. 



This is to avoid the maximal value of 9 which leads to vanishing dshi ft ■ 
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For all the cases we find dja ^ 10~^. The superluminal motion is realized when dose > 
\dmass\ and this condition can be satisfied for all existing setups. However, as follows from 
the Table I, oscillations can not explain the OPERA result: Indeed, we have shown that 

• the additional distance of the i/^— propagation is too small: It is restricted by size of 
the neutrino wave packet. Since the signal in OPERA is not suppressed, P/^/^ ~ 1, 
the effect can not be related to oscillations into sterile neutrinos. Therefore for known 
Am^ the distance is further suppressed by the phase 0p acquired by neutrinos in the 
production region. 

• the distance has strong dependence on neutrino energy E; d increases as E decreases. 

In MINOS at certain energies increase of propagation distance can be realized due to 
smallness of P^^. In the case of Idecay -C ^ 1) the effect can be enhanced if one 

includes in consideration the exponentially suppressed tails of the wave packets utilizing Ip 
as the wave packet formation region instead of Idecay (Of course, this requires an extremely 
intense beam of neutrinos.) From eq. f l55]) we obtain 

So, in comparison to the case y = 1, for y ^ 1 the distance increases by a factor 
r)s{oo)/rjs{l) = 12.71. Therefore in all the cases in the Table, where y = 1 the distance 
would be at most 12.71 times larger. In this way we obtain dose ~ 0.02 cm for OPERA 
(1 GeV), MINOS (1 GeV) and T2K (0.4 GeV) instead of numbers in the Table I. Even 
with such an enhancement we have dosc/c < 10~^, and observation of these distances looks 
practically impossible. 



B. Enhancement of contributions 

Let us consider a possibility to increase the additional distance of z/^ propagation and 
to observe effective superluminal motion of z/^. With nanotechnology the accuracy of time 
measurement would allow probing changes of distances of propagation as small as 1 m. 
Notice that one can evaluate the change of time of i/^ propagation for a given baseline L 
dividing the corrections to the distance of propagation by v, e.g. Atosc = —dose/v ~ —dose/c. 

In this discussion we will use the results for box-like shape factor which are simple and 
transparent, and at the same time are valid upto 0(1) coefficients also for the exponential 
shape factors. The result for the box-like wave packet fHTl) can be rewritten in the following 
form 

/ \ / , xSin^2^ . / , 0„\ . ^, 

{x)osc = -f^^(0p)^T75 — sm + — , (54) 



2P V 2 



where 



2 



sin ^ — ^ "^os ) , (55) 



andP^^ = P^^(0 + f: 



According to (|5^ there are several ways to increase dose ^ {x)c 
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1). Increase a = a{T,ldecay,v.„,6^^). For this one should lower the neutrino energy, 
use neutrino parents with large lifetimes, increase the size of decay region and use off-axis 
neutrino beams. According to |15[ 



where 9^^ is the angle between a neutrino momentum in the pion rest frame and a momentum 
of pion in the laboratory frame. For 9^^ = it reproduces our previous results. With increase 
of 9^^ the size of the wave packet increases. Simultaneously, the energy of neutrino decreases 
in such a way that axE is invariant [l5|. E.g. for 9^^ = 7r/2, ax = and energy becomes 2 
times smaller. Thus, for off-axis experiments the additional distance of propagation can be 
larger. For the backward emitted neutrinos, 9^^ = tt, the size of the wave packet becomes 
the largest one: ax = 27^f r° = 2t>r^. Then ax = 2lp, if Idecay > Ip- The energy becomes 
times smaller, thus for E^^ = 1.4 GeV it will be about 2 MeV. 

2) . Increase of the phase (pp = lp{lp,ldecayJu) acquired over the neutrino production 
region. This is possible with increase of Am^, decrease of the neutrino energy, etc.. 

3) . Selection of certain value of the oscillation phase, l,^), by the selecting particular 
values of the baseline or /and neutrino energy. Varying the oscillation phase we find that 
maximal value of the distance equals to 

{x{t))ZT = -Sa g{<Pp) , =, (57) 



2^{ci + s^y-[sm^29'-^^y 
where S = sign[sm{(l) + (pp/^)] and it is achieved at determined from the equation 



, , ^„ , sin^ 29 / sin ^ 
cosU+^ =-— — — ^ . (58) 



2 J c^ + s^ 

4). Select baseline and/or neutrino energy so that P^^ is small. 

Here we consider two possibilities to increase dose ~ {x(t))osc- (i) Decrease of the survival 
probability P^^ which implies strong suppression signal, (ii) Increase of a, g{(f)p) and sin0 
without substantial decrease of P^^ and therefore the signal. Let us explore them in order. 

1). Small P^^ can be achieved by selecting certain values of (j){L,E). We consider the 



/ifi 

case of small (pp in which 



Notice that according to f HOj) the minimum of P^^ is obtained at = vr — 0p/2 and equals 

P™" = cos2 20 + sin2 29-^. 

fj,fj. 

However, for this value of the distance dose vanishes to the second order in (pp. On the 
other hand, for (p = tc we obtain 

P;";" = cos^2^ + sin2 2^^, 
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and the distance equals 



dn 



sin^ 29 0-0- 
cos2 2^ + sin^ 48 



(59) 



In the case of maximal mixing this equation gives dose = cr/4. 
In the limit (^p — )■ we obtain from (|57l) 



imax 
^osc 



crc 



which corresponds to 



COS( 



For maximal mixing we find from ([60 



sin^ 26 



W2^ + §sin^2^ 



(60) 



2s^c 



2„2 



d 



max 
osc 



a 



2v^' 



which can be considered as the maximal possible additional distance of propagation due to 
oscillations. This, however, corresponds to very small survival probability: P„„ = -^f sin^ 29. 



For large 0p and = vr we obtain 



sin^ 29 



2cos2 2^ + 0.5sin2 2^ (^1 



sin (pp 
<j>p 



COS I 



sm( 



and therefore dn 



a/2. 



2). Let us maximize other factors in 
reaches maximum 



MD for P^^ = 0(1) . The function g{ 



0.2 



at (hn determined from the condition 



tan(0p/2) 



bj2 



1 - 



The smallest value of the angle which satisfies this equation is 



TT + e, where e > 0. 



(Other values give local extrema.) With further increase of (j)p (above vr) the function g 
decreases as l/0p, and for small (pp it has linear dependence: g oc 0p. 

The oscillatory factor in f l5^ gives maximal distance of propagation when (p + 0p/2 = 
37r/2 so that sin(0 + 0p/2) = —1. This, in turn, can be achieved by selecting the baseline 
or/and neutrino energy. For (pp ^ tt which corresponds to maximal value of g, one needs 
to have = vr + 2TTk {k is an integer). For these phases the probability equals the average 
probability: P^^ 



+ s^. Consequently, in the case of maximal mixing we obtain 



ag{'K) 



0.2a. 



(61) 



For fixed value of (f)p the contribution {x{t))osc can be slightly larger if + 0p/2 deviates 
from 37r/2. For (pp = n and maximal mixing the equation (|57|1 gives {x{t))osc = 0.26(T. At 
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the same time the probabihty becomes smaller than before: P^^ ~ 0.3. Essentially we have 
here enhancement of the additional distance due to decrease of probability. 

Thus, the only way to get significant additional distance of i/^ propagation (which equals 
significant fraction of the length of the wave packet size, {x(t))osc/<^) without strong sup- 
pression of signal is to increase the phase acquired by neutrinos over the production region, 
so that (j)p = 0(1). 

Let of consider possible setups which can realize 0p = 0(1) and increase a. Apparently 
the region of formation of neutrino wave packet is restricted by the length of decay tunnel. 
At most /p = (1 — 2) km and therefore according to fl33|) we have two possibilities: 

1) . For known Am^ = 2.5 ■ 10~^ eV^ one can take low energies: E = few MeV. For 
muon neutrinos this can be arranged with slowly moving muons which have two orders of 
magnitude larger lifetime. The detection should be then via the neutral current interactions. 
Alternatively, one can detect backward moving neutrinos from the pion decay. 

For the electron neutrinos one can use not only muon decays but also nuclei decays (beta 
beams). I'e detection via the charged currents interactions is of course easier. In this case, 
however the effect is suppressed by small 1-3 mixing: sin^ 2^13 < 0.1. This allows one to 
achieve {x{t))osc/<^ ~ (0.01 — 0.1). At low energies the size of the wave packet can be large: 
cr ~ (0.1 — l)/p, and therefore, {x{t)) 

Qsc ~ (1 — 10) m. However, realizations of these setups 
will still require significant loss of signal (selection of low energy part of the muon decay 
spectrum, use the backward neutrinos, small cross-sections, etc.). 

2) . If sterile neutrinos exist with large mass splittings, Am^ ~ 1 eV^, the phase 0p can 
be large at higher energies: ~ 1 GeV. Here again the effect will be suppressed by small 
allowed mixing: sin^ 26^14 < 0.1. As a result, {x(t))osc/<^ ~ (0.01 — 0.02) can be achieved. 

For low energy pions the size of the neutrino wave packet and the neutrino formation 
region are determined by the decay length, and therefore 

Am2 Am'^E^ g 
0. = -^^--^-^ = (62) 

This phase weakly depends on energy and turns out to be of order unity at Am^ > 10 eV^. 
The length of the wave packet equals a = r°/27^ and for E,, = 0.5 GeV we obtain a ~ 1 
m. This gives {x{t))osc ~ (1 ~ 2) cm. 

Larger distance can be achieved for neutrinos from the muon decay. In the GeV energy 
range Idecay > Ip- So, according to (1221) a ~ (10 — 20) m, and {x(t))osc ~ (10 — 20) cm. 

Summarizing, the large additional distance of z/^ propagation, dose ~ (0.1 — 10) m, can be 
obtained in rather non-standard experimental setups with low energy accelerators of muons 
and long decay tunnels, with use of the electron neutrinos or in the presence of large Am^. 



VI. OSCILLATION IN MATTER AND NEUTRINO VELOCITIES 

In most of experiments, including OPERA, neutrinos propagate in matter. In this con- 
nection let us discuss influence of the matter effect on the neutrino velocity [13]. Notice 
that in the absence of mixing the dispersion relation in matter reads 

E = + ^2 + 

where V is the matter potential which can be written at low energies as Vq ~ /Sy^Gpn with 
Gp and n being the Fermi constant and the electron number density in matter, respectively. 
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The constant /3 depends on the neutrino flavor. Since at low energies the potential Vq does 
not depend on momentum of neutrino, dVo/dp = 0, the group velocity remains unchanged 
in matter: 

dE _ p 
'dp ~ E' 

V depends on energy due to the W boson propagator. For the elastic scattering in forward 
direction = 0, and therefore the energy dependence of V appears when the W exchange 
occurs in the s-channel. In usual media this is possible for Uf, only when z/g annihilates with 
electrons. In this case 

2 2 



TTiyy — s — 2meP — ml ' 
and now 

dV _ 2mem^ 

dp (m^-s)2' 

If s < m'^ we obtain the energy independent contribution to velocity: 

A.„«!!!!i!^L = l,8xlO-»f^)f-^), (63) 
'< V3g/cmV VO-5/ 

where p is the matter density and Y^. is the electron fraction. This conclusion does not 
depend on whether V is constant function of x or not. The contribution to the neutrino 
velocity in (l63l) is too small to affect any of our discussions, and hence will be ignored. 

Since neutrinos are mixed, the propagating degrees of freedom are neutrino eigenstates 
in matter. The energy of these states (eigenstates of the Hamiltonian) are given by 



^u=P+^ + ^±5\/(''-2^-2») +4(— ) sln^2e. 

Differentiating by p we obtain 

^1,2 ~ A^l I /l„2 I 7, \^^) 



4p2 4p2 



(cos2^^-|^) +sin2 2^ 



In the limit of zero potential or small energies it reproduces the usual result. At the resonance 
point, where the denominator is minimal, we find 



«M-l-^T^s„>2.. (65) 



At very high energies or large matter potential the velocities equal 

^ m? + mi Am^ , , 

<2-l-^4^±^cos2^^. (66) 

Thus, as follows from flMj) the correction to the velocity due to mixing in matter alone 
cannot lead to superluminal motion. 
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For the — Ur mixings in the hmit of zero 1- 3 mixing (neglecting loop corrections) 
the difference of potentials is zero, = 0, and the situation is reduced to the vacuum 
case described in Sec. II - IV. With non-zero ^13 in the three neutrino case the eigenvalues 
acquire additional dependence on momentum related to the matter potentials. Furthermore, 
the mixing angle and energy splitting in matter should be taken into account. However, at 
E > 6 GeV, i.e. above resonance energy in the Earth, this dependence is weak, and it can 
be neglected in the first approximation. 

It is possible to extend these statements onto the case of three flavor neutrino mixing. 



The energy eigenvalues in matter can be written as E" 
given by Is!' 



p + 



2p 



1, 2, 3), where Aj are 



Ai 
A2 
A3 



1 



3t 



u 



s vs^ 



3t 



u 



± r 
-s + -uy t 



3t, 



(67) 



with 



s - 
t -- 

u 



A21 + A31 + a, 
A2iA3i + a[A2i(l • 
'1 _i /2s3 



S?2Cy +A3l(l 



COS 



-COS 



9st + 27aA2iA3icf2C?3 

2(s2 - 3t)3/2 



(68) 



In ( 168|) we have used the notations Ajj 
of Aj is only through a. Therefore, 



m| and a 



Vp. Notice that p-dependence 



A,; 



2p^ 



^ 1 dXi da 
2p da dp 



+ 



1 dXi 



2p2 



d(\og Xj 
d{\oga) 



(69) 



According to flB7|) (this feature can be seen in the plot for Aj as function of a given, e.g. 



m 



19l . |20|) Aj is a monotonically increasing function of a, but with the growth rate slower 
than a. Therefore, ^jj^^J]^^ < 1- Asymptotically, the largest eigenvalue A3 behaves as a, 
approaching to the equality. Thus, f j — 1 < which excludes the possibility of neutrino's 
superluminal velocity due to matter effect. 



VII. CONCLUSION 

There are several factors which change the shape of the wave packet of a muon neutrino, 
and consequently, influence the distance of propagation for a given time, and hence the 
velocity of neutrino: (i) the relative shift of the wave packets of the mass eigenstates, (ii) 
oscillations, (iii) absorption, (iv) production. In this paper, we focussed on the first two 
factors which are mutually correlated: both are due to the mass squared difference, and 
therefore the separation of the wave packets is always accompanied by oscillations and vice 
versa. 

1. We have computed the distances of z/^ propagation in the presence of mixing and 
oscillations. The oscillations lead to distortion of the shape factor of the z/^ wave packet. 
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This, in turn, changes the effective distance traveled by neutrinos and therefore the group 
velocity. This is essentially related to the oscillation effect within the neutrino wave packet. 
The oscillatory pattern is squeezed and therefore the effect is enhanced in the same way as 
the size of the neutrino wave packet shrinks in comparison to the pion decay length or size 
of decay tunnel. 

We find that the distance of the z/^ propagation is proportional to the length of the 
wave packet and the oscillation phase 0^ acquired by neutrinos along the decay path of 
the parent particles (pions, K-mesons ) where neutrino wave packet is formed dose oc (70p. 
Furthermore, dose has an oscillatory behavior with distance determined by the oscillation 
length. For small distances, L < 1^/2, the oscillations reduce the group velocity. The 
distance dose becomes positive for baselines L = l„/2 — etc.. In this range of baselines 
motion can be effectively superluminal. The additional distance is restricted by the size of 
neutrino wave packet: dose < (^xu- The additional distance strongly decreases with increase 
of energy: d^sc ~ 1/-^^. 

2. Distortion of the z/,,— wave packet is also produced due to relative shift of the wave 
packets of the mass eigenstates even in the case when oscillation effect does not depend 
on coordinate. This effect previously considered in the literature is proportional to Avt — 
Aw? /2EH, and therefore neghgible in comparison with the oscillation effect. 

3. For the OPERA setup with Am^ = 2.5 ■ 10"^ eV^ and ~ 17 GeV we obtain 
dose ~ —10^^ cm and other contributions are vanishingly small. Therefore the OPERA 
result, which corresponds to the distance ~ +20 m, can not be explained. For average energy 
E — YJ GeV the oscillations reduce the distance propagated by neutrinos and consequently 
the group velocity. Furthermore the distance dose strongly decreases with increase of the 
neutrino energy. 

4. We estimated the additional distances of propagation for different experimental 
setups. In particular, we find dose = (0.01 - 0.04) cm for MINOS and T2K. Change of the 
time of propagation (for a given L) can be obtained as At ~ —dose/c- 

5. Larger additional distance can be obtained for neutrinos from decays of particles 
with longer lifetime: muons, nuclei (beta beams). The dose becomes larger with decrease of 
neutrino energy and for large Am?, if exist. It can be as large as several meters. However, 
this requires rather non-standards experimental setups and huge neutrino fluxes. 

6. Measurements of additional distances of the flavor neutrino propagation opens a way 
to determine sizes of neutrino wave packets and test certain quantum mechanical features 
of neutrino production. In any case, the oscillation effect should be taken into account in 
analysis of future measurements of the neutrino velocities. 
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Appendix A: Maximal effect of shift of the eigenstates 



Here we clarify the effect of shift of the wave packets discussed in the hterature. FoUowing 
01, Isl, 01 we take Ap = 0, so that the oscillation phase is the same along the wave packet 
{i.e. by itself oscillations do not produce distortion of the packet). For simplicity we take 
the box-like shape factors for the mass eigenstates. 

Let us consider evolution of neutrino state produced as z/^. According to ([T]) in the state 
the z/i and z/2 packets have amplitudes c and s respectively. (Here we omit normalization 
factors which are irrelevant for final result.) The amplitude of probability to find z/^ in the 
ui part of the state is and in the 1/2 part: s^. Result of measurement of z/^ in the full 
evolved neutrino state is then determined by interference of the z/^ parts of ui and z/2. The 
interference pattern, in turn, is determined by the relative (oscillation) phase of the two 
mass eigenstates, and is described by cos The phase does not depend on x being the 
same over whole wave packet. 

At the moment of time t after production the centers of the z/i and 1^2 packets have 
coordinates x = vit + a/2 and x = V2t + a/2 correspondingly. The average value of the 
two coordinates \s x = vt = (t>i + V2)t/2 + a /2. The relative shift of the packets equals 
Aft = /\m?t/2E'^ . Correspondingly, the front edges of vi and V2 are at vit + a and V2t + a. 

Consider the shape of the wave packet and iV'i^pP- As we marked in Sec. IVC, due 
to the relative shift of the wave packets of mass eigenstates there are three different spatial 
parts of the wave packet {yi > V2): 

1. In the front edge part, {v2t + a) < x < {vit + a) only Ui packet is present, so that 
ipi^^ = and l'?/'jy^P = c^. The size of this part is given by the shift Aft. 

2. In the overlapping part, (fit) < x < (f2t + a), both wave packets are nonzero and 
li^iyj"^ = + COS$p. 

3. In the trailing edge part, (f2t) < x < (fit), only z/2 wave packet is present and 
= s^, SO that \ipuj'^ = s^. 

In the first approximation the distance of z/^ propagation is determined by the position of 
the "center of mass" of the wave packet squared. In the case of maximal mixing whole the 
picture is completely symmetric with respect to x and therefore (x) = x. If mixing deviates 
from maximal one the shift of the packets leads to asymmetric distortion of the z/^ wave 
packet. This in turn, leads to a shift of the "center of mass" from x: 

(x) = X + 6x. 

For non-maximal mixing (with c > s) the forward edge is higher by c*^ — = cos 26 than 
the trailing edge. This difference should be compensated by the shift of the center of mass 
by amount Sx in the overlapping region. This compensation leads to the condition 

Aft cos 261 = 26x\c^ + s^cos0|l 

Therefore 

_ cos 26* Aft 

dx 



2\c^ + cos 0p 

The maximal shift would correspond to cos*!* = —1. It implies the destructive interference 
in the overlapping region, when = tt, so that l^pu^l"^ = — = cos^ 29, and therefore 

cos 26* , 1 , . , 

Sx = Aft- — - = Aft- -, Al 

2 COS"' 2(J 2 cos 26 
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which reproduces results in j^, [6|, [7|. Thus, a superluminal motion here is a result of 
interplay of the coordinate-independent oscillations and the relative shift of wave packets 
due to different group velocities. Apparently the shift 5x is restricted by the size of the wave 
packet. In this case it is clear that v has no physical meaning. It is the velocity of the 
"center of mass": some effective point in the fiat overlapping part of the shape factor. No 
one real structure in the shape factor (edges of different regions) is moving with v > c. 
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